Draft version August 1 1 , 2008 

Preprint typeset using I4TgX style emulateapj v. 10/09/06 



THE KINEMATICS OF KEPLER'S SUPERNOVA REMNANT AS REVEALED BY CHANDRA 

Jacco Vink 1 

Draft version August 1 1, 2008 

ABSTRACT 

I have determined the expansion of the supernova remnant of SN1604 (Kepler's supernova) based on archival 
Chandra ACIS-S observations made in 2000 and 2006. The measurements were done in several distinct energy 
bands, and were made for the remnant as a whole, and for six individual sectors. The average expansion 
parameter indicates that the remnant expands on average as r oc t 5 , but there are significant differences in 
different parts of the remnant: the bright northwestern part expands as r oc i° 35 , whereas the rest of the 
remnant's expansion shows an expansion r oc t 6 . The latter is consistent with an explosion in which the 
outer part of the ejecta has a negative power law slope for density (p oc v~ n ) of n = 7, or with an exponential 
density profile^ oc exp(— v/v e )). The expansion parameter in the southern region, in conjunction with the 
shock radius, indicate a rather low value (< 5 x 10 50 erg) for the explosion energy of SN1604 for a distance of 
4 kpc. An higher explosion energy is consistent with the results, if the distance is larger. 

The filament in the eastern part of the remnant, which is dominated by X-ray synchrotron radiation seems to 
mark a region with a fast shock speed r oc t 7 , corresponding to a shock velocity of v — 4200 kms -1 , for a 
distance to SN1604 of 4 kpc. This is consistent with the idea that X-ray synchrotron emission requires shock 
velocities in excess of ~ 2000 km s -1 . 

The X-ray based expansion measurements reported are consistent with results based on optical and radio 
measurements, but disagree with previous X-ray measurements based on ROSAT and Einstein observations. 

Subject headings: X-rays: observations individual (Kepler/SN 1604) -supernova remnants 



1. INTRODUCTION 

Among the more than 250 known Galactic supernova rem- 
nants, the remnants of the historica l supernovae hold a special 
place (IStephenson & Gree n 2002). This has partially to do 
with the fascination for historical events that caught the imag- 
ination of astronomers in ancient China, the Middle East, and 
in renaissance Europe. But also more scientific reasons make 
the study of historical supernova remnants worthwhile: we 
know the exact age of the objects. Moreover, the historical 
supernova remnants are among the youngest supernova rem- 
nants, which means that their X-ray emission is largely dom- 
inated by shock heated ejecta rather than shocked interstellar 
matter. Historical remnants are therefore prone to offer new 
insights into the supernova explosion properties. 

The y oungest historical superno va remnant (SNR) is 
SN1604 (Steph enson & Greenl |2002|) . also known Kepler's 
SNR (Kepler for short). 2 The supernova was first sighted 
on the evening of October 9, 1604, low above the horizon. 
It owes its early discovery probably to the much anticipated 
simultaneous conjunction of Jupiter, Saturn and Mars. Jo- 
hannes Kepler lived in Prague at that time, and suffered from 
bad weather. However, from the first reports on he took a keen 
interest in the new star, and after the weather improved, he 
started his own observations. The results of his observations 
and his correspondence with other observers led to his book 
on the supernova, "De Stella Nova", which was published in 
1606. 

The SNR of SN 1604 has been a puzzling object for some 
time (Blair 2005). Both the historical light curve of the su- 
pernova "(iBaaddl 119431) and its relatively high Galactic latitude 
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2 Cas A is a younger remnant, but not st rictly an historical remna nt, be- 
cause the supernova was likely not observed I Stephenson & Green 2002). 



(I = 4.5°, b = 6.4°) suggest that Kepler is the result of 
a Type la supernova. However, optical observations of the 
remnant reveal the presence of copious amounts of nitrogen, 
in particular in the Northwest. Nitrogen is an element as- 
sociated with s tellar winds rath er than Type la supernovae. 
This prompted Bandiera (1987) to suggest that the progeni- 
tor was a massive, runaway star, thus explaining both the ori- 
gin of nitrogen and the high Galactic latitude of the super- 
nova. Nevertheless, the X-ray spectrum of Kepler indicates a 
large abundance of iron, which points to a Type la supern ova 
dKinugasa & Tsunemil [19991: ICassam-Chenai et al.ll2004 . A 
recent deep Chandra observation revealed no evidence for a 
neutron star dReynolds et alJl20O7l) . which would be expected 
if SN1604 w as a core collaps e super nova. Having considered 
all evidence, iReynolds et al.l (120071) conclude that Kepler is 
the remnant of a Type la supernova from a relatively massive 
progenitor star, attributing the nitrogen to stellar wind loss 
from either the white dwarf progenitor, or the companion star. 

The situation c oncerning the ki nematics of Kepler is 
equally confusing. Dick el et all (119881) studied the expansion 
of Kepler in the radio, based on VLA 6 cm and 20 cm data 
covering a time span of 4 yr. They found expansion rates in- 
dicating significant deceleration with, on average, an expan- 
sion rate consistent with a radial expansion law r oc t°, 
and in the northern part even as low as r oc f 35 . Op- 
tical expansion measurements, based on ground based and 
Hubble Space Telescope imaging of the bright Ha filaments 
in the Northwest covering a ti me span of 16.33 y r, showed 
proper motions of 1.3" to 1.6" (ISankrit et alJl2005l) . This cor- 
responds to an expansion law following r oc t° 35 , consis- 
tent with the radio measurements. However, X-ray expansion 
measurements, based on Einstein and ROSAT observations 
with the high resolution imagers onboard the se two satellite s 
indicated nearly free expansion: r oc < 93 (Hughes 1999). 
Note that of all wavelength regimes the X-ray emission is 
most closely associated with the dynamics of the remnant, 
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TABLE 1 
Observations 



Observation 


Start date 


MJD 


Exposure (ks) 


116 


2000-06-30T12:04:56 


51725.5034 


49.45 


6714 


2006-04-27T23: 13:49 


53852.97 


159.84 


6715 


2006-08-03T16:52:09 


53950.70 


161.17 


6716 


2006-05-05T19:18:18 


53860.80 


160.05 


6717 


2006-07-13T19:04:25 


53929.79 


108.18 


6718 


2006-07-21T14:55:14 


53937.62 


109.18 


7366 


2006-07-16T05:24:07 


53932.23 


52.12 



since the shock heated plasma has most of its emission in 
X-rays, and best compares with hydrodynamic simulations. 
The optical emission is confined to a region close to the shock 
front, whereas the kinematics as derived from the radio emis- 
sion should be related to the X-ray emitting plasma, but the 
mismatch between radio an d X-ray expansion measurements 
for Kepler, but also Cas A dVink et al.ll 19981: iKoraleskv et all 
1998[lDelaney & Rudnickll2003l) . and Tycho's SNR dHughesI 
20001) complicates the interpretation. 

Resolving the discrepancy between radio and X-ray ex- 
pansion measurements, and obtaining the overall kinemat- 
ics of Kepler's SNR is important for several reasons: First 
of all the distance to the SNR is poorly known. One way 
to estimate the distance is to measure both the proper mo- 
tion and the shock velocity. The latter can be done inde- 
pendent of the proper motion by measuring the broad Ha 
emission line, which is the result of charge transfer from 
neutral hydrogen entering the shock to shock heated pro- 
tons. The broadening of the line therefore is a measure of 
the proton temperature, which is linked to the shock veloc- 
ity. For Kepl er this method yields a shock veloc ity of 1550- 
2000 kms- ^Fesen et alJll989tlBlair et al.lll99"ll) . implying a 
distance of 3.9tj g kpc dSankrit et al.ll2005l) . This is slightly 
lower, but consistent with the most recent me asurement based 
on H I absorption features, 4.8 ± 1.4 kpc dReynoso & Gossl 
fl999l) . 

Another issue is the dynamical state of the SNR. As long as 
the mass of the ejecta dominates over the the mass of shock 
heated circumstellar material the remnant is little decelerated, 
and said to be in the free expansio n phase with R oc t (see, 
however iTruelove & McKed ll999). Once the energy in the 
shock heated circumstellar medium dominates the total en- 
ergy, the SNR is said to have entered the Sedov-Taylor stage 
of its evolution with R oc t aA for a uniform density medium. 
If SN1604 is a Type la supernova one expects it to have a rel- 
atively low ejecta mass of 1.4 AT©. A high expansion rate, as 
found in X-rays, is therefore puzzling. 

Finally, in recent years it has been found that all young 
SNRs have thin filaments at the shock front, whose emission 
is dominated by synchrotron radiation (e.g. IVink & Lamingl 
l200l IBamba et al.1 l200l IVolk et al.1 l2005lh The widths 
of these filaments can be used to infe r magnetic fields 
dVink & Laming! 120031: IBamba et al.ll2005t IVolk etail 120051: 
IWarren et al.ll2005T e.g.), which turn out to be relatively high. 
This is therefore e vidence for cosmic r ay driven m agnetic 
field amplification dBell & Lucekl l200lt iBelH 12004b . How- 
ever, it is not quite clear how the magnetic field scales with 
density and shock velocity; it cou ld be eithe r B 2 oc pv 2 
(IVolk et alJl2005h or B 2 oc pv 3 JBel 12001: Ivirikl l2006h . 
The range in densities among the SNRs is quite large (a fac- 
tor of 100 from SN1006 to Cas A), so the dependency of 




Energy (keV) 

FIG. 1. — The Chandra A CIS-S X-ray s pectrum of the entire remnant, 
based on all 2006 observations I Reynolds et al. 2007). The arrows and dotted 
lines mark the energy bands used to measure the expansion. 

B 2 on p can be determined reasonably well. However, X- 
ray synchrotron fila ments only arise for high s hock velocities 
v > 2000 kms~ 1 (Ahar onian & Atovaniri999h . and since all 
known SNRs have v < 6000 kms -1 , the dynamic range in 
velocity is not very high, as compared to the dynamic range in 
densities. Moreover, the uncertainties in the measured veloc- 
ities is quite high. A more accurate assessment of the shock 
velocities in thos e regions where X-ray synchrotron filaments 
have been found (Reynolds et al. 2007) is therefore valuable. 
In addition, the cut-off photon energy of X-ray synchrotron 
radiation depends n ot only on v s , but also on th e cosmic ray 
diffusion parameter (iAharonian & Atoyanll 19991) : 

with ?/ a parameter that relates the actual diffusion coeffi- 
cient to the diffusion coefficient for the most optimal diffu- 
sion coefficient for fast cosmic ray acceleration, the so called 
Bohm diffusion. For Bohm dif fusion r\ = 1, and the magnetic 
field is highly turbulent (see Malk ov & Drurv||200lL for a re- 
view). Apparentl y, r? « 1 for young SNRs (jyin^|200l|2Q06|; 
IStage et al]|2006l) . Accurate shock velocity measurement are 
therefore important for estimating the diffusion constant, and 
the related turbulence of the magnetic field. 

Here I present expansion measurements based on archival 
Chandra data obtained in 2000 and 2006. 

2. DATA AND METHOD 
2.1. Observations 

Chandra observed Kepler several times, the first time in 
June 2000 as part of the Guaranteed Time Observation pro- 
gram. The most rece nt observation was ma de in 2006 as part 
of the Large Program dRevnolds et al.l2007l) . All observations 
were made with the ACIS-S array, with the SNR being pro- 
jected on the ACIS-S3 chip. 

Since proper motion measurements are more reliable for 
longer time spans I limited the analysis to the 2000 and 
2006 observations, ignoring an additional observation made 
in 2004. The 2006 observation was split in several pointings, 
which we list in Table Q] Weighing the Modified Julian Dates 
(MJD) of the 2006 with the exposure times I obtain an average 
MJD for the 2006 observations of 53912.13. I therefore find 
an average time span between the 2006 and the 2000 Chandra 
observations of At = 5.985 yr. 
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For the analysis I used the standard processed ("evt2") 
events files obtained from the Chandra data archive. All event 
files were processed by the Chandra data center in 2007. 

2.2. Method 

The expansion measurements were made using the same 
method and updated C ++ code employed for the SNR Ca s 
A using Einstein and ROSAT observations (IVinket al.ll 1998b . 
The results for Cas A were later confirmed by measurement s 
based on Chandra observations dDelanev & Rudnickl 120031) . 
The method is similar to what has b een used to m easure the 
previous X- ray expansion of Kepler dHughesll 19991) . and Ty- 
cho's SNR (Hughes 2000). One aspect not explored by these 
studies is the dependence of the expansion on the energy band, 
which was not possible due to lack of energy resolution of 
Einstein and ROSAT. 

For the current analysis I extracted images in several en- 
ergy bands using custom-made software, which, as described 
below, allows one to center the image on a given sky co- 
ordinate and correct for bore sight errors. I chose energy 
bands based on the spectroscopic features of the X-ray spec- 
trum (see Fig. Q}: 0.5-0.7 keV (covering the O VII/O VIII 
line emission), 0.7-1.0 keV (Fe XVII-Fe XXI L-shell emis- 
sion), 1.0-1.5 keV (Fe XXII-Fe XXIV L-shell emission, per- 
haps blended with Ne IX/X and Mg XI/XII emission), 1.7-1.9 
keV (Si XIII K-shell emission), 2.0-4.0 keV (covering K-shell 
emission from Si, S, Ar, and Ca), and 4.0-6.0 keV, which is 
dominated by continuum emission (a combination of thermal 
bremsstrahlung and synchrotron radiation). 

The expansion factor, / was measured by regridding the 
images from the 2006 observations, which are statistically su- 
perior, using a simple expansion law: 

x 2 = round{ x + f(x\ - x n ) + Ax} (2) 
t/2 = round{ y Q + f(y 1 - y Q ) + Ay}. 

with (xi,yi) being the original pixel coordinates, (3:2,2/2) 
the new pixel coordinates (rounded to the nearest integers), 
and (Ax, Ay), free parameters which are fitted to correct for 
pointing errors between the observations. The absolute roll 
angle accuracy is very high ~ 25", and therefore role an- 
gle errors do not contribute to registration errors. 3 Note that 
pointing errors are indistinguishable from errors in the ex- 
pansion center. For the expansion center (xo, 2/0) I adopted 
an estimate of the geometrical center of the remnant based 
on the X-ray image: a/2000 = 17 h 30 m 41.25 s and <5j2ooo = 
-21°29'32.95". 

For obtaining the best fit expansion rates, maxi mum likeli - 
hood statistic for Poisson distributions was used (ICashll 19791) . 
which minimizes: 
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with riij being the counts in pixel of the ObsID 116 
(2000) image, and rriij being the predicted counts based on 
the 2006 image, which has been regridded (Eq. [3j and renor- 
malized such that E^n^- — 'Sijiriij. The statistical fluc- 



tuations are dominated by the observation in June 2000, with 
its 49 ks exposure. I can therefore treat the combined 2006 
observations with an exposure of 750 ks as a model. Note 
that the last term in Eq. [3]can be ignored for fitting purposes, 
because it depends only on riij which does not change from 

3 See http://cxc.harvard.edu/cal/ASPECT/roll_accuracy.html, and Tom 
Aldcroft (private communication). 
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FIG. 2. — The expansion parameter as function of energy. The overall 
expansion rates, based on fitting the expansion for the whole remnant, are 
indicated in black. The errors correspond to 90% confidence ranges. Due 
to the brightness of the northwestern region, the expansion rates are biased 
toward the expansion in the Northwest. The expansion rate averaged over 
all 6 sectors are displayed in red, with errors corresponding to the standard 
deviation. 



one set of fit parameters to the other, as long as the numbers 
of bins over which is summed remains unchanged. The er- 
ror in the fit parameters can be estimated using the fact that 
AC = C — C mm is similar to A x 2 (ICashll 19791) . The fitting 
procedure itself is done by scanning the relevant parameter 
range, iteratively switching between the various parameters, 
and with each iteration decreasing the step size. 

The code has the option to fit only certain regions of the 
image, using a combination of image masks and region files. 
For the overall fits a mask was used, based on the broad band 
image, blocking out all pixels falling below a certain thresh- 
old level. I fine tuned this mask, such that the SNR image, 
smoothed with a = 2 pixels plus a border of 2 additional 
pixels, falls within the mask. 

A problem with Eq. [3]is that one can only sum over model 
pixels rriij that are not zero. Regridding the model image 
makes that in low emissivity regions a model pixel may acci- 
dentally be set to zero for one set of parameters, in which case 
it is ignored, whereas for other values of the fit parameters it is 
non-zero. Hence the number of pixels over which the statis- 
tic is derived is not constant. This problem does not occur 
for Kepler's SNR for energy bands with sufficient statistics, 
such as those covering the bright Fe-L emission, but it is im- 
portant for the continuum image. In order to overcome this 
problem the model image (based on the 2006 observations) 
was smoothed with a Gaussian with a — 1 pixel, ensuring 
that the number of pixels over which the statistic is derived 
remained constant. I checked for the images with the best 
statistics, whether smoothing had any effect on the expansion 
measurements, but within the statistical error the small scale 
smoothing did not affect the measured expansion rates. For 
that reason I adopted Gaussian smoothing to all 2006 images, 
in order to have one consistent way of measuring for all en- 
ergy bands. 

The Chandra ACIS chips have a pixel resolution of 0.492", 
slightly undersampling the telescope point spread function. 
Because of bore sight effects absolute coordinates are accu- 
rate up to ~ 0.4". This means that by adding all the 2006 
observations one may introduce a slight blurring by approxi- 
mately 1 pixel. In order to start from the best possible images 
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I used the same code as for expansion measurements reported 
below, but fixing / = 1, and using a broad band (0.3-7 keV) 
image. All fits were made with respect to the ObsID 6714 im- 
age. Having fitted (Ax, Ay) for each individual observation, 
the final extraction of the images was made with corrections 
for the individual bore sights, after which images in the same 
energy band were added together. The average bore sight was 
0.4 pixels in both coordinates, with rms errors of 0.2 pixels. 

Throughout this paper I use three different ways of char- 
acterizing the expansion of Kepler's SNR: 1) the expansion 
rate defined as R = (1 - f)/At, with At = 5.985 yr; 
2) the expansion time r exp = l/R, which is perhaps the 
most intuitive number, as it gives the age of the remnant in 
case one assumes free expansion; 3) the expansion parameter 
P = T SNi604/Tea:pj the ratio of the true age of Kepler's SNR 
over the expansion time. 

From an hydrodynamical point of view (3 is the most 
important parameter. In general the shock radius, r s , of 
SNRs in distinct differen t phases evolves with time r as 
dTruelove & McKedl 19991) : 



TABLE 2 

Expansion rates 



with K a constant. This gives for the shock velocity: 



T 



this shows that 
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with d the distance, and the angular radius. 

3. RESULTS 

3.1. The average expansion 

The simplest approach to measuring the expansion of Ke- 
pler is to fit for each energy band the expansion factor and 
bore sight/expansion center offsets (Ax, Ay). Of course both 
the expansion center and bore sight errors should not depend 
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FIG. 3. — Confidence ranges for the expansion rates for the remnant as 
a who le, based on the log likelihood ratio, AC, with respect to the best fit 
(Cash 1979). From left to right the curves correspond to the energy bands 
0.5-0.7 keV (O VII/O VIII line emission), 0.7-1.0 keV (FeL line emission), 
1.0-1.5 keV (Fe-L line emission), 1.7-1.9 keV (Si XIII line emission), and 
4.0 - 6.0 keV (continuum emission). 





R 




Ay 


o / A / A . r. J\ 

ri (Axy i\y fixed) 




(vo/yr) 


pixels 


pixels 


(wyr) 


0.5-0.7 keV 


0.084 ± 0.004 


-0.58 


-0.38 


0.084 ± 0.004 


0.7-1.0 keV 


0.108 ±0.002 


-0.58 


-0.38 


0.106 ±0.002 


1.0-1.5 keV 


0.115 ±0.002 


-0.67 


-0.45 


0.108 ±0.002 


1.7-2.0 keV 


0.119 ±0.003 


-0.57 


-0.38 


0.117 ±0.003 


2.0-4.0 keV 


0.131 ±0.004 


-0.74 


-0.47 


0.129 ±0.004 


4.0-6.0 keV 


0.158 ±0.006 


-0.70 


-0.61 


0.155 ±0.007 


Mean (± rms) 


0.119 ±0.025 


-0.64 


-0.44 


0.116 ±0.024 




0.48 ±0.10 






0.47 ±0.10 




840 ± 181 






859 ± 185 



on the energy band. So a second iteration involves fixating 
(Ax, Ay), and then measuring the expansion of the remnant. 

However, still this is not ideal; Kepler's SNR is bright in 
the northwest, so expansion measurements of the whole rem- 
nant are skewed toward the northwest. Moreover, since a next 
step involves measuring the expansion in different regions, er- 
rors in (Ax, Ay) results in errors in the expansion rate, with 
opposite signs for opposite sides of the SNR. 

I therefore adapted a different method: I first divided the 
SNR in 4 sectors of each 90°. Instead of fitting the four 
sectors individually, I paired them in North-South and East- 
West pairs. For the North-South pair I fixated Ay and for 
the East- West pair Ax. This ensures that there is no inter- 
ference between proper motion and bore sight corrections, 
since (Ax, Ay) are measured perpendicular to the proper mo- 
tions. This procedure was repeated twice in an iterative way 
for the four energy bands with sufficient counts (the oxy- 
gen, 2 Fe-L bands, and Si band). The best fit values are 
(Ax, Ay) = (-0.48 ± 0.01, -0.35 ± 0.04), which were then 
used for all subsequent expansion measurements. 

Table [3] lists the expansion rates for the individual energy 
bands both with (Ax, Ay) as free parameters, and fixed to 
the aforementioned best fit values. Fig. [2] shows the expan- 
sion parameter j3 for the individual energy bands. Comparing 
the expansion rate measurements for fixed and fitted Ax/ Ay 
shows that the measured expansion rates are, within the sta- 
tistical error, identical. The likelihood ratios (AC) for the 
individual fits show smooth curves as a function of expansion 
rate, indicating that the minimum value for the statistic is well 
defined, with no sub-minima (Fig. [3]). 

There is a clear tendency for the expansion parameter to in- 
crease with the photon energy , ranging from (3 = 0.34 ± 0.02 
for the oxygen band, to /3 = 0.62 ± 0.03 for the continuum 
band. Note that these values are significantly lo wer than the 
expansion parameters reported by Hughes] (Il999h . 

Fig.[4]shows the effect of taking the average expansion into 
account for the 1.0-1.5 keV band, by showing the difference 
between the 2000 and 2006 images. 

The difference image shows some image artifacts accross 
the northern shell. These are caused by streaks in the June 
2000 images. The streaks are not immediately obvious in 
the images themselves, but the difference image emphasizes 
them, in particular when the expansion has been corrected for. 
Similar artifacts s howed up in measuring the expansion of Cas 
A with Chandra (iDelaney & Rudnickll2003h . From the resid- 
uals after expansion correction, I estimate that the brightness 
errors accross the streaks are typically ~ 10%, except near the 
bright knot, where the brightness error peaks to ~ 30%. Here 



FIG. 4. — The difference between the 2000 and 2006 1.0-1.5 keV image. A small gaussian smoothing has been applied to the images (<r = 1 pixel). The image 
on the right has only been corrected for the bore sight error. The 2006 image on the left has also been corrected for the av erage expansion. Note that the left 
image still shows some expansion of the filament in the east, and overcorrects some of the expansion in the west (see section 1X2) . Also note the image artifacts 
in the north, which are more visible in the left hand image, because there the expansion has been taken out, leaving the artifacts as one of the dominant sources 
of differences between the two images. 



is not quite clear whether the large deviations are caused by 
the streaks or by the knots itself. I come back to the streaks 
below. 

3.2. The expansion as a function of azimuth 

Given the strong asymmetry of Kepler's SNR it is quite nat- 
ural to expect that th e expansion will al so be asymmetric, as 
indeed was found by Dic kel et al.l {1988) in the Radio. 

I investigated this by dividing the remnant in six sectors, 
each spanning 60°in azimuthal angle (Fig.[5]l. In order to mea- 
sure only proper motions around the shell of the remnant, the 
central region was ignored, which features a bar-like structure 
of unknown origin (Fig. [5]). For the measurements of the ex- 



FIG. 5.— Color image combining the 1.0-1.5 keV (red), 1.7-2.0 keV 
(green) and 4-6 keV (blue) with overlayed the sectors used for measuring 
the expansion as a function of azimuth. The circle indicates the region not 
considered for the expansion measurements. 



pansion per sector, the same expansion center and bore sight 
corrections were applied for each sector, as for the mean ex- 
pansion measurements reported in section [3~T1 

In total 36 expansion rates were measured, covering six sec- 
tors and six X-ray bands. The expansion rates as a function of 
azimuthal angle are listed in Table [3] whereas the expansion 
parameters and expansion times are displayed in Fig. [6] The 
expansion rates averaged over all sectors show less variation 
as a function of energy than the mean expansion rates based 
on the whole remnant. An important difference is that the 
averaged rates give equal weight to all sectors of the remnant, 
whereas the mean expansion rate is biased toward the brighter 
northwest of the remnant (Fig. O. 

Table [5J and Fig. [6] show that the expansion in the north- 
northwestern sectors is considerably slower than in the other 
parts of the remnant: for the northwestern sector the expan- 
sion parameter ranges between j3 — 0.3 and (3 = 0.4, corre- 
sponding to T exp = 1000 — 1500 yr. For the southern and 
eastern sectors this is (3 = 0.55 — 0.68, corresponding to 
Texp — 590 — 730 yr. Fig. [6] also suggests that the varia- 
tion in the average expansion as a function of energy (Fig. [6j» 
can be mostly attributed to the northwestern region. 

As shown in Fig. [4] there are streaks in the June 2000 im- 
ages, which happen to be confined to the northern sector. The 
expansion measurements in this sector gives an expansion pa- 
rameter similar to that of the neigboring northwestern sector. 
Given the presence of artefacts in this region one should treat 
this expansion measurement with more caution than those of 
the other 5 sectors. In order to have a quantitative estimate 
of the streaks on the expansion measurements of the north- 
ern region, I also measured the expansion rate after blocking 
out most of the pixels affected by the streaks, thus removing 
about half of the northern sector. Using this smaller region 
resulted in a higher expansion rate, corresponding to j3 = 0.5, 
except for the 4-6 keV band, which remained at (3 = 0.3. 
This difference could be either attributed to the removal of the 
streaks, i.e. the value listed in Table|2]for the northwestern re- 
gion could be affected by systematic errors of order 30%, or 
it could be that the regions blocked out have a slower expan- 



6 



TABLE 3 
Expansion rates per sector 







0.5-0.7 keV 


0.7-1.0 keV 


1.0-1.5 keV 


1.7-1.9 keV 


2.0-4.0 keV 


4.0 - 6.0 keV 


Sector 




R (%/yr) 


R (%/yr) 


R (%/yr) 


R (%/yr) 


R (%/yr) 


R (%/yr) 





(N) 


0.074 ± 0.007 


0.098 ± 0.002 


0.096 ± 0.003 


0.105 ± 0.005 


0.108 ± 0.006 


0.080 ± 0.015 


60 


(NE) 


0.094 ± 0.011 


0.133 ± 0.004 


0.141 ± 0.005 


0.134 ± 0.007 


0.140 ± 0.010 


0.171 ± 0.021 


120 


(SE) 


0.168 ± 0.022 


0.164 ± 0.010 


0.174 ± 0.006 


0.160 ± 0.011 


0.173 ± 0.008 


0.178 ± 0.014 


180 


(S) 


0.155 ± 0.028 


0.148 ± 0.012 


0.164 ± 0.013 


0.143 ± 0.014 


0.154 ± 0.014 


0.109 ± 0.016 


240 


(SW) 


0.116 ± 0.018 


0.160 ± 0.008 


0.143 ± 0.009 


0.139 ± 0.011 


0.143 ± 0.013 


0.173 ± 0.021 


300 


(NW) 


0.068 ± 0.006 


0.077 ± 0.003 


0.081 ± 0.004 


0.093 ± 0.007 


0.103 ± 0.009 


0.138 ± 0.021 


mean (± rms) 


0.113 ± 0.042 


0.130 ± 0.035 


0.133 ± 0.037 


0.129 ± 0.025 


0.137 ± 0.027 


0.142 ± 0.04 



0.8 - 



u 

<D 
*J 

a 

CO 
CO 

a 

c 
o 

n 
C 
CO 

a 
x 



0.6 - 



0.4 



0.2 




100 200 
Position Angle 



likely that most, if not all, of the continuum emission from 
this region is thermal bremsstrahlung. The most unambigu- 
ous X-ray synchrotron emitting filament is the arc-like fila- 
me nt in the east. As show n by [Cassam- Chenai et all {2004) 
and iRevnolds et alj d2007l) . the spectrum of this filament is 
completely dominated by continuum e missio n. The width of 
this filament was used bv IVolk eT al. (2005) to infer a mag- 
netic field of 250 fiG (see also Bamba et al. 2005). The 2006 



FIG. 6. — The expansion parameter of Kepler's SNR as a function of azimuthal angle. The angle is measured from the north in a counter clockwise direction. 
The different colors indicate expansions measured in different energy bands, using the same color coding as in Fig. [5] The data points have been cyclic, so the 
data points at 360° repeat those at 0°. The vertical axis on the right indicates the corresponding expansion time. 

sion rates. The latter option is quite well possible, since the 
streaks affect mostly the western part of the northern sector, 
i.e. the region closests to the more slowly expanding north- 
western sector. In that case the different expansion are due to 
a gradient accross this sector, which is quite plausible. Given 
the uncertainty, I will not explicitely discuss the northern re- 
gion, but instead concentrate on the contrast between the slow 
expansion in the northwest and the rest of the remnant. 

3.3. The X-ray synchrotron filament in the Southeast 

As mentioned in the introduction, the shock velocity near 
X-ray synchrotron filaments is of considerable interest for 
both understanding magnetic field amplification (BelJ 120041 : 
IVolketalJl200H IVinkN2006l) and the magnetic field turbu- 
lence. 

Like Cas A and Tycho, Kepler's SNR shows contin- 
uum emission around the whole perifery of the remnant, 
some of which is probably X-ray synchrotron emission (e.g. 
ICassam-Chenai' et al.ll2004l) . However, in the northwest the 
continuum emission is more diffuse and associated with re- 
gions of the most intense line radiation. It is therefore quite 



ever, n ote that the estimate of the magnetic field bv lVolk et alj 
(120051) assumes spherical symmetry in order to estimate the 
true, deprojected, filament width. Taking the 2.5"at face value 
results in a mag netic field of 85 fiG using the formula in 
IVink et aTl (120061) . This should be taken as a conservative 
lower limit. 

Since it is important to know under what conditions X-ray 
synchrotron emission occurs, I measured the expansion of the 
eastern filament separately. The spectrum from this fil ament 
shows hardly any line emission (IRevnolds et al J 120071) . so I 
measured the expansion using a 0.3-7 keV broad band image 
in order to improve the statistics of the measurements. The 
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FIG. 7. — Left: Detail of the eastern part of Kepler's SNR, showing in red the mask used for determining the proper motion of the filament, in green the broad 
band image, and in blue a smoothed version of 4-6 keV image. Right: Emissivity profile of the northern region of the filament, based on the 2006 Chandra 
observations in the 4-6 keV band. 

expansion rate that is found is R = (0.176 ± 0.007)%, cor- 
responding to a relatively large expansion parameter of (3 = 
0.71 ± 0.03. The filament is located at an angular radius of 
2.1' from the center. Since X-ray synchrotron filaments trace 
the shock front, one can therefore translate the expansion pa- 
rameter into a shock velocity of v s = (4200±170)<i4kms _1 , 
with d,4 the distance in units of 4 kpc. This is twice as fast as 
the shock velocities inferred from optical sp ectral and proper 
motions studies in the northwestern region dBlair et all 1991 
Sankrit et al. 2005), but consistent with the value adopted by 
VolketalJ (1207)1 7 

Since the shape of the filament has a radius of curvature 
smaller than the radius of the remnant, one may wonder in 
what directions the filament is actually expanding: Is the fil- 
ament caused by a blow out, in which case one expects the 
expansion center to be closer to the approximate curvature 
center of the filament, or is the expansion center close to the 
geometrical center of the whole remnant? 

In order to get some handle on this, I also fitted the expan- 
sion, but leaving the center of expansion as free parameters. In 
that case, the best fit center of expansion was more toward the 
west (ajsooo = 17 h 30 m 40.77 s , (5j 20 oo = -21°29'29.49") 
than the adopted center, i.e. opposite of the non-thermal fil- 
ament. This, and in addition blinking of the 2000 and 2006 
images by eye, suggests that the curved structure is moving 
more or less as a coherent structure, rather than expanding 
from a center close to the filament. This is reminiscent of the 
kinematics of a bow shock structure. The expansion param- 
eter did not change substantially, when the expansion center 
was treated as a free parameter: = 0.67 ± 0.04. 

4. DISCUSSION 

I have measured the expansion of Kepler's SNR us- 
ing archival Chandra data from observations performed in 
2000 and 2006. These new X-ray expansion measurements 
largely agree with expansion measurements based on radio 
dDickeletal.il 19881) and optical (ISankrit et al.ll2005l) measure- 
ments. Specifically, the results confirm that the average ex- 



pansion parameter is (3 « 0.5. The expansion as a func- 
tion of azimuthal angle shows a clear difference in expansion 
rate between the northwestern and other parts of the remnant 
with the northwestern part having an e xpansion parameter 
(3 w 0.3 - 4, as also found in the radio (iDickel etalJI 19881) 
and optical (Sankrit et al, l2005l) . and the other parts having 
j3 » 0.6, in agreement with the radio measurements. 

The expansion measurements presented here are in dis- 
agreement with previous X-ray measur ements, based on 
ROSAT and Einstein data dHugheslll999l) . which suggested 
(3 « 0.9. In terms of resolution this new measurement should 
be better than the ROSAT-Einstein measurement, despite the 
long baseline of the latter, 17.5 yr. The resolution of the high 
resolution imagers on board the ROSAT and Einstein satel- 
lites is about 4", amounting to 0.22" yr -1 for a baseline of 
17.5 yr. The pixel resolution of the Chandra ACIS instrument 
is 0.43", so the resolution per unit time for the present study 
is about 0.07" yr" 1 . It is difficult to assess what causes the 
discrepancy between the new re sult and the ROS AT/Einstein 
study, as the measurements by Hughes] d 19991) used a simi- 
lar method as em ployed here, and for the expansion of Cas A 
dVinket al.ll 1998L also based on Einstein and ROSAT data). 

4.1. The expansion parameter in theoretical models 

Theoretical models of SNR evolution predict d irectly the 
expan s ion parameter, a dimen s ionles s number (jChe valier 
1982; iDwarkadas & Chevalierl 119981: iTruelove & McKee 
19991) . The best known example is the evolution of the shock 
wave in the so called Sedov-Taylor phase, which treats the 
supernova as a point explosion in a uniform density medium. 
This gives (3 = 0.4. IChevaherl d 1982b analyzed the early evo- 
lution of a SNR, in the context of power law ejecta density 
models, i.e. p oc v~ n . This gives j3 = (n — 3) /n for an explo- 
sion in an uniform density medium. For Type la supernovae 
it has been argued that n = 7, which should therefore result 
in (3 = 0.57 during the early phase of the SNR. 

The expansion parameter in this case refers to both the con- 
tact discontinuity (separating shocked circumstellar medium 
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(CSM) and shocked ejecta) and the forward shock. The 
plasma in the shell has a range of expansion parameters cen- 
tered around (3 — 0.57, ranging from B — 0.63 near the re- 
v erse shock, to (3 = 0.43 n e ar the forward shock. 

iDwar kadas & ChevalieJ (Il998h studied the hydrodynami- 
cal evolution of Type la remnants numerically, using both 
power law ejecta density profiles, and exponential ejecta den- 
sity profile (p oc exp(— v/v e )). They explicitly provide the 
expansion parameter of the shock itself, which in case of 
the exponential density profile ranges from (3 — 0.8 very 
early in the evolution to j3 = 0.4, although never reaching 
(3 = 0.4. The expansion parameter of the shocked plasma 
can be obtained from their Fig. 3: At a late phase, around 

t « 500Ef^ 5 n Q 1 ' 3 yr (uq being the preshocked number 
density and E$i the explosion energy in units of 10 51 erg) for 
both types of density profiles (3 « 0.5 — 0.6 for the forward 
shock, but for the plasma the expansion parameter is almost 
uniformly (3 w 0.38. The lower expansion parameter for the 
plasma is not surprising, given that the plasma directly behind 
the shock moves with v — 3v s /4 in the case of a monatomic 
gas with adiabatic index 7 = 5/3. Therefore, in the late phase 

ftshock — 4/3p/ asma /3. 

These studies are important for interpreting the measured 
expansion parameters of Kepler's SNR. First, it is good to be 
aware that the expansion parameters of the shock may be dif- 
ferent from the plasma. In the late phase, the plasma expan- 
sion parameter is lower than that of the shock itself, whereas 
in the early phase of the evolution there is a range of values, 
but around the contact discontinuity the expansion parameter 
is similar to that of the shock. Therefore, during the early evo- 
lution the expansion parameter is expected to be close to the 
expansion parameter of the forward shock. 

4.2. Inferences on the shock velocities 

The question now arises what the X-ray expansion mea- 
surements really provides: the expansion parameter of the 
shock or that of the plasma behind it? The shock velocity 
is a pattern speed, and this is certainly part of what is mea- 
sured. As Fig. |4] shows, correcting for the mean expansion 
removes the strong fringes in the difference map, around the 
forward shock. On the other hand, also the velocity of the 
plasma itself may influence the best fit expansion parameters. 
The measurements themselves are skewed toward the outer 
part of the shell, simply because the proper motions are larger 
there, and more pixels are involved. 

The models discussed above therefore suggest that the mea- 
sured expansion parameter is a lower limit to the expansion 
parameter of the forward shock, late r in the evolution of the 
SNR (i.e. t' > 1 in the notation of Dwark adas & Chevalier ! 
fT998L their Eg. 6). In an earlier phase, for which the lChevalierl 
U982J) models may be applicable, the measured expansion pa- 
rameter may be a good representation of the expansion param- 
eter of the forward shock. 

The expansion measurement for most of the remnant, ex- 
cept for the northwestern part, gives (3 « 0.6, which is con- 
sistent with the expansion parameter of {3 = 0.57 for both 
the contac t discontinuity a nd the forward shock for the n = 7 
model of I Chevalier! (119821) . the prefered model for Type la 
supernovae. It is also in the approximate range for the for- 
ward shock expansion param eter for the exponential density 
profile model considered by Dwarkadas & Chevalier! (119981) . 
but only in case the remnant is still in the early phase of its 
evolution, i.e. if the self-similar time variable is t' < 1. 



Note that the theoretical values for the expansion parame- 
ter, whether for the shock or the plasma, are all lower than 
the expansion parameters obtained from the previous X-ray 
measurements based on ROSAT and Einst ein observations 
(!Hugheslll999h . namely 8 = 0.93. The IChevalien (fl982h 
model for n — 7 gives for the plasma velocities a range of 
(3 = 0.4 — 0.6, whereas the exponential density profiles gives 
(3 < 0.8 after the first few de cades in the life of the SNR 
(IDwarkadas & Chevalierlfl998l) . In addition,the previous X- 
ray expansion measurements are inconsistent with the Fe-K 
emission detected to come fro m reverse shocked ma t erial i n 
Kepler's SNR, as discussed by ICassam-Chenai' et al.1 (f2004): 
Fe-K emission requires a well advanced reverse shock and 
high reverse shock speed, which is incompatible with free ex- 
pansion. 

The value for the expansion parameter in the southern and 
eastern sectors are consistent with the e xpans ion parame- 
ter for the n — 7 model of Chevalier (1982), which ap- 
plies to both the contact discontuity and the forward shock. 
It seems, therefore, justified to translate the measured ex- 
pansion parameter into a shock velocity, using that the av- 
erage angular radius of Kepler is r = 1.76'. This gives 
v s = 3010(/3/0.6)e?4 kms -1 , with e?4 the distance in units 
of 4 kpc. 

For the northwestern region shock velocities have been 
measured using the widths of Ha lines, which gives a 
direct measurement of the post-shock proton temperature. 
These suggest shock velocities in this region of 1500- 
2800 kms- ^Fesen etalJIigsgt lBlairetal.ll 19911). with the 
most likely value centered around 1660 km s _1 (ISankrit et al.l 
120051) . The pro per motions of the H a emitting shock regions 
is 0.088" yr -1 (ISankrit et al.ll2005l) . This corresponds to an 
optical expansion parameter of (3 = 0.33. This value is con- 
sistent with the value reported here for the X-ray expansion 
parameter for the northwestern region, (3 — 0.3 — 0.4, and 
with the radio expansion parameter j3 « 0.35. 

Finally, the highest expansion reported in this paper is for 
the X-ray synchrotron filament in the east {(3 — 0.71). The 
size of the X-ray synchrotron emitting region is determined 
by the loss time of the highest e nergy electrons and their ad- 
vectio n away of the shock front ( Vink & Laming 20031 IVinkl 
2006). As a result the size of the X-ray emitting region is 
expected to be fixed as a long as the magnetic field and the 
shock speed are approximately constant. The displacement 
of the filament therefore reflects the speed of the shock front, 
rather than the movement of the plasma. This suggests that 
the measured filament movement is a pattern movement, and 
reflects directly the movement of the shock front, implying 
v s = 4200d 4 kms" 1 . 

This relatively high shock velocity and the protrusion of 
the filament outside the general shock radius of Kepler's SNR 
make it quite likely that the filament marks a part of the shock 
that is exp anding into a low dens ity region, as previously sug- 
gested by ICassam-Chenai' et al.l (120041) . Moreover, the high 
velocity and the low density both help to explain why in this 
region X-ray synchrotron emission dominates over thermal 
emission. The high velocity is consistent with the idea that 
only shocks with velocities v s > 200 km s" 1 give rise to 
detec t able X-ray synchrotron radiation (Aharonia n & Atovanl 
[T999l:IVink et al]l2006t iHelder & Vinkll2008l) ~ 

4.3. The expansion and the energetics ofSN1604 

The regions of Kepler's SNR that are likely to be most re- 
vealing about the explosion parameters of SN1604 are the 
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FIG. 8. — The deprojected 1-1.5 keV (Fe-L) image of Kepler's supernova 
remnant, used to estimate the location of the reverse shock (indicated by the 
dotted line with a radius of 1.3'). The deprojection was made in 18 inde - 
pendent sectors, following the procedure described in Helder & Vink 12008). 
The brightness scale is in percentage per bin, scaled in such a way that the 
total adds up to 100% integrated flux in each sector. (Figure kindly provided 
by Eveline Helder.) 

southern/southwestern regions. In the northwestern region it 
is clear that t he expansion para mete r is too low to fit either 
the mo dels of IChevaiien (1 19821) and iDwarkadas & Chevalier] 
( 1998), or the Sedov evolution. In the east the radius of the 
remnant is not well determined, due to the protruding non- 
thermal filament. For the southwestern region the shock ra- 
dius is 1.77', corresponding to a physical radius of r s = 
2.06^4 pc. As an additional constraint we can use the ra- 
dius of the reverse shock. From the 1-1.5 keV (Fe-L) im- 
age the typical reverse shock radius is estimated to be ~ 1.3' 
(Fig. [HJ, corresponding to a reverse shock radius of r r w 
1.6^4 pc. Using these values, and the measured expansion 
parameters for the southern/south western regions (0 ~ 0.6 
in conjunction with the models of iTruelove & McKee (1 19991 
or IDw arkadas & Chevaherl (119981) . we can constrain the pre 
shock density and explosio n energy for Kepler' s SNR . 

Fig. [9] shows the n = 7 [Truelove & McKed (1 19991) model 
for two different choices of kinetic energy and circumstellar 
medium density, assuming a distance of 4 kpc. One choice 
is to assume that rig. = 1 cm~ 3 , the other that the energy 
is the canonical explosion energy of Eq = 10 51 erg (= 1 
Bethe). In both cases it has been assumed that the ejecta mass 
is M e j = 1.4 Mp,, as is to be expected for a T ype la su- 
pernova (IWooslev et al.ll2007t iMazzali et al]|2007l for recent 
discussions on this issue). It is clear that both in terms of the 
reverse shock position in Kepler's SNR and in terms of mea- 
sured expansion parameters a lower kinetic energy is to be 
prefered. 

Comparing the measured properties for the southern 
part of Kepler's SNR with the numerical models of 
IDwarkadas & Chevalier! (1 19981) leads to a similar conclu- 
sion. The measured expansion parameters are consistent 
with a low value of the sel f-similar time coordinate in 
IDwarkadas & Chevalier! d!998l) . i.e. t' < 1, with (3 w 0.6, 
corresponding to t' « 0.6. However, a conservative up- 



per limit is t' < 2, corresponding to f3 = 0.5 for the 
shock itself. Translating the dimensionless time, and the as- 
sociated normalized shock radii (0.85 < r' < 1.6) into a 
physic al ages and shock radius using the conversion equa- 
tions in Dwarkadas & Chevalier ( 1998, Eq.4-6), one finds that 
0-7 < n H /(cm- 3 ) < 5, and 0.2 < E /(1B) < 0.5, with a 

preference for the lower v alues. 

So for both the ITruelove & McKed d!999l) and 
IDwarkadas & Chevalier! (119981) models it appears that 
SN1604 had a relatively low explosion energy. How- 
ever, this depends also on the adopted distance of 4 kpc 
(ISankrit et a . 2005). For 5 kp c , close to the nominal distance 



estimate of iReynoso & Gossl (1 19991) . the allowed range of 
densities and energies is 0.4 < fin/fcm" 3 ) < 2.5, and 
0-3 < E /(1B) < 0.9, which is still rather low compared 
to the uniform kinetic energies of 1.2 ± 0.2 B i nferre d 
for observed Type la supernovae (Wooslev et al. 2007). 
Moreover, the upper limit requires a relatively high density 
for a SNR located 470gLi pc above the Galactic plane. Only 
for an adopted distance of 6 kpc is the angular radius of 
the forward and reverse shock of the SNR consistent with a 
kinetic energy of 10 51 erg. The associated ISM density is 
then Tin = 0.5(cm~ 3 ). A distance considerably further than 
4 kpc seems therefore preferable, provided that SN1604 was 
indeed a Type la supernova, and considering the evidence that 
most Type la superno vae have energies in excess of 10 51 erg 
(IWoosley et al.ll2007l) . 

A similar conclusion was recently obtained based on the 
non-detect ion of TeV emission fr om Kepler with the H.E.S.S. 
telescope (A haronia n et al.l2008l) . Note that the conclusion of 
lAharonian et all d2008l) depends on assumptions concerning 
the TeV luminosity of the remnant, which depends on the ex- 
plosion energy and on the fraction of the energy that goes into 
accelerating cosmic rays. 

4.4. Estimates of the swept up mass in the Northwest 

The northwestern region of Kepler's SNR has an expansion 
parameter even lower than the expansion parameter expected 
for the Sedov-Tay lor phase. This is probably caused by a non- 
uniform density profile. The northwestern part of the remnant 
shows in the optical nitrogen-rich material, suggesting a shell 
ejected by the progenitor system. The fact that the expansion 
rate in the northwestern region seems to be a function of pho- 
ton energy may be accounted for if the inner layer of shocked 
ejecta is hotter, and if in this layer the expansion parameter 
is larger. The slowest expansion is measured for the oxygen 
band, and it has been suggested that most oxygen emission 
is fro m shocked circumstellar med ium, rather than from the 
ejecta (Cassam-Che nai et al.l l2~004). The trend with energy in 
this region could be an effect of a slow response of the in- 
ner layers to density enhancements closer to the shock. This 
reinforces the idea of a dense shell that has recently been en- 
countered by t he blastwave. 

As noted by Reyn olds et al.l (120071) the presence of nitro- 
gen rich material from the progenitor system is difficult, but 
not impossible, to reconcile with SN1604 being a Type la su- 
pernova. A possibility is, for example, tha t the supernova 
belon ged to the short-lived Type la channel ( Mannucci et al. 
2006), in which the white dwarf progenitor and its companion 
star, were relatively massive stars. However, for the remnant 
of SN1604 the complication remains that the progenitor, or 
its companion, must have deposited a substantial amount of 
mass at a large distance. 

The amount of material must have been substantial in order 
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FIG. 9. — The evolution of shock radii (left) and expansion parameters of the foward shock (right) according to the n = 7 model of Truelove & McKee 
1 1999) for two different sets of kinetic energies and circumstellar densities. The values were adapted such that the forward shock matches 2.06 pc (indicated 
by the upper dotted line) at the age of Kepler's SNR (402 yr), valid for a distance of 4 kpc. In the top panels riu = 1 cm" 3 , resulting in a kinetic energy of 
Eq = 2 X 10 50 erg. In the bottom panels Eq = 10 51 erg, in which case nu = 12 cm -3 . In this case both the position of the reverse shock (the lower dotted 
line in the left panels) and the expansion parameter do not match the measurements. 

for the shock to have decelerated so much that the expansion 
parameter is even lower than expected for the Sedov-Taylor 
phase. If we assume the material encountered in the north- 
west was a shell, covering a fraction of about / = 0.25 one 
can estimate the mass in the shell by requiring that the mass in 
the shell must be more than the swept up interstellar medium 
in other parts of the remnant, which have nn ~ 1 cm" 3 . This 
corresponds to a mass of M swept = l(//0.25)nH<i4 M Q . It 
therefore seems reasonable to assume that the mass encoun- 
tered in the northwest is also about 1 M Q . 

This material must have been lost from either the progenitor 
of the supernova, or from its companion star. This suggests 
a strongly non-conservative binary evolution scenario: Too 
much mass in the shell means less mass available for accre- 
tion onto the white dwarf, complicating its evolution toward 
a Type la supernova. Another problem may be how to eject 
this material to a distance of of ~ 2 pc from the progenitor. 
Perhaps, the shell is caused by nova explosions on the progen- 
itor? 

The remnant of SN1604 remains, therefore, a puzzling, but 
intriguing object. Due its unusual properties it may in the fu- 
ture reveal new aspects of Type la supernovae. Future studies 
of the SNR may provide some answers, but a high priority 
would be to identify light echos of SN1604, and obtain their 
optical spectra. This has b een done very recently for the SNRs 
Cas A (K rause et al.l2 008), which appears to have been a Type 
lib supernova, and SNR 0509-67.5, an energetic Type la su- 



pernova in the Large Magellanic Cloud (Rest et al. 2008). 

5. SUMMARY 

The expansion of Kepler's SNR (SN1604) was measured 
using archival Chandra data. The expansion in all parts of the 
remnant is inconsistent with free expansion with the expan- 
sion measurem ent of (3 ss 0.9, that was previously reported 
(lHugheslll999l) . 

The X-ray measurements reported here, and previous radio 
and optical expansion measurements, show that the remnant 
expands more slowly in the bright northwestern part, (3 rs 
0.3—0.4, than in the rest of the supernova remnant, where (3 « 
0.6. The fastest expansion is found for the X-ray synchrotron 
filament in the eastern part of the remnant, (3 « 0.7. 

The remnant seems not yet to have entered the Sedov- 
Taylor phase of its evolution: Apart from the northwest- 
ern region the expansion parameters are consisten t with the 
early e x pansion phase as de t ailed in the models bvlChevaliei 
1982); [Truelove & McKed d!999l) : iDwarkadas & Chevaliei 
1998). For the northwestern part of the remnant a differ- 



ent scenario is needed, since it's expansion parameter is even 
smaller than for a Sedov-Taylor evolution model. 

The kinematics of Kepler's SNR, in particular the more 
undisturbed southwestern region, is only consistent with cur- 
rent hydrodynamical models of Type la SNRs, if its distance 
is cons iderably large r than the 4 kpc o b tained bv lSankrit et alj 
( 120051) . Recently, lAharonian et all ([2008) suggested that 
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d > 6 kpc, based on the non-detection of Kepler by the 
H. E.S.S. TeV 7-ray telesc ope. Both in the present paper and 
in lAharonian et all (|2008), the conclusion regarding a large 
distance is based on the assumption that Kepler was indeed a 
Type la SNR with an explosion energy > 10 51 erg. 

The X-ray synchrotron filament seems to move with 
a sho ck speed of 4200^ 4 kms -1 , consistent with the- 
orv (lAharonian & Atovanl 1 19991) and other observations 
(Held er&Vinkl 120081) that indicate that only shocks with 



> 2000 km s 1 emit X-ray synchrotron radiation. 
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